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Strange Electric Form Factor of the Proton 
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By combining the constraints of charge symmetry with new chiral extrapolation techniques and 
recent low-mass quenched lattice QCD simulations of the individual quark contributions to the 
electric charge radii of the baryon octet, we obtain an accurate determination of the strange electric 
charge radius of the proton. While this analysis provides a value for G S E {Q 2 = O.lGeV 2 ) in agreement 
with the best current data, the theoretical error is comparable with that expected from future 
HAPPEx results from JLab. Together with the earlier determination of G S M , this result considerably 
constrains the role of hidden flavor in the structure of the nucleon. 



One of the great challenges of modern hadron physics 
is to unravel the precise role of hidden flavors in the struc- 
ture of the nucleon. Because of their relatively light mass, 
strange quarks are expected to play the biggest role and 
it is with respect to strangeness that there has recently 
been enormous experimental progress. Indeed, new re- 
sults on strangeness in the nucleon have been reported 
recently from the HAPPEx and GO Collaborations 
at JLab, which complement earlier work at MIT-Bates 
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and Mainz The situation is by far the best at Q 2 = 
0.1 GeV 2 , where in addition to the usual linear combina- 
tion of electric and magnetic form factors, a measurement 
of parity violation on 4 He allowed an accurate extraction 
of G S E , namely G S E {Q 2 = 0.1 GeV 2 ) = -0.013 ± 0.028. 
A new experimental investigation, in which this error is 
expected to be reduced by roughly a factor of two, was 
conducted in late 2005. This makes it imperative to find, 
as we report, a way to make a theoretical estimate of 
comparable accuracy — or better. 

With respect to the strange magnetic from factor of 
the proton we recently reported a calculation an order of 
magnitude more precise than current experiments 0, . 
This calculation exploited the advances in lattice QCD 
which have enabled quenched QCD (QQCD) simulations 
of magnetic moments at pion masses as low as 0.3-0.4 
GeV 0, HI @, as well as the development of new chi- 
ral extrapolation techniques 0, [nj . To minimize the- 
oretical uncertainty, an essential input was the precise 
(experimental) measurements of the magnetic moments 
of the ground state hyperons — a luxury unfortunately 
not available for charge radii. Nevertheless, we show here 
that even the limited data on hyperon charge radii, when 
combined with new lattice simulations and chiral extrap- 
olation techniques, yield a precision commensurate with 
the published data. Alternatively, using the best esti- 
mates for the charge radii of the valence u and d quarks, 
extracted from the lattice simulations, yields a determi- 
nation of the strange quark contribution to the proton 
charge radius with an uncertainty comparable with that 
anticipated from the latest HAPPEx measurement. 
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FIG. 1: Diagrams illustrating the two topologically different 
insertions of the current within the framework of lattice QCD. 



As illustrated in Fig. 2] the three point function re- 
quired to extract an electromagnetic form factor in lattice 
QCD involves two topologically distinct processes (each 
incorporating an arbitrary number of gluons and quark 
loops). The left-hand diagram illustrates the connected 
insertion of the current to one of the "valence" quarks of 
the baryon. In the right-hand diagram the external field 
couples to a quark loop. The latter process, where the 
loop involves an s-quark, is entirely responsible for G S E . 

Charge symmetry, which is the invariance of the QCD 
Hamiltonian under rotations by 180 degrees about the 
2-axis in isospace, is typically respected at a level better 
than 1% in nuclear physics [12j. This gives one confidence 
in using it to extract the strange form factors experimen- 
tally and we shall also use it. Under charge symmetry, 
the electric or magnetic form factors of the octet baryons 
satisfy |l3j: 
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e u u p + e d dP + O n . 
e d u p + e u dP + O n . 
e u u s + e s s s + O s 



(1) 
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Here, p, S _ etc, represent any form factor of the physical 
proton and X - , and similarly for the other baryons. For 
the present application we take them to be the mean- 
square charge radii, defined as — 6 d/dQ 2 Ge{Q 2 )\q^=o- 
Similarly, vP denotes the contribution to the proton 
mean-square charge radius from the two valence w-quarks 
(for u-quarks of unit charge), see LHS of Fig.Q] 
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Of course, the strange quark contribution is entirely 
contained in the quark-loops, On, illustrated on the right 
hand side of Fig. However, one has to separate the 
contributions associated with u and d quarks from those 
involving s quarks. To this end, we define u, l d and 
l s as the mean square charge radii of the loop contribu- 
tions associated with u, d and s quarks (of unit charge), 
respectively. Hence we may write: 
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(2) 



where the ratio of s- and d-quark loops, R s d = l s/ l d, is 
expected to lie in the range (0,1). Note that, in deriving 
Eq. (J2J, we have used charge symmetry (or m u = m,j) 
to set e u = e d 0| . Since the chiral coefficients for the d 
and s loops in the RHS of Fig. 1 are identical, the main 
difference comes from the mass of the K compared with 
that of the ir. 

With a little algebra On, and hence e s, may be isolated 
from Eqs. (JTJ and J5J: 



2p + n - -= (E+ - E" 



(3) 



It is worth emphasizing that this expression is an exact 
consequence of QCD, under the assumption of charge 
symmetry. 

There is a relation similar to Eq. J3J involving the 
charge radii of the 3 baryons and the ratio of the charge 
radii of u-quarks in the n and 3 , but as there is no ex- 
perimental data at all for the cascades we do not show 
it. Even in the case of the E hyperons, there is data only 
for the E~ — from the SELEX Collaboration at Fermi- 
lab Contrasting the magnetic moments, the experi- 
mental error is relatively large, E~ = — 0.61±0.12±0.09 
fm 2 . Nevertheless, as we shall see, even this is enough, in 
combination with lattice QCD estimates of the ratio of 
charge radii in the E + compared with the E _ , to obtain 
an accurate estimate of the strangeness radius. 

Combining the experimentally measured mean square 
charge radii of the p and n (0.757 ± 0.014 fm 2 and 
—0.116 ± 0.002 fm 2 , respectively [l5j) we find: 
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0.61 ±0.12 ±0.09 



(4) 



To obtain On, the necessary input from lattice QCD is 
now the two ratios in Eq.Q, namely and E + /E~. 

Figures |2 El and |5] show the state of the art calcula- 
tions from which we can extract these ratios. The nu- 
merical simulations of the electromagnetic form factors 
presented here are carried out using the Fat Link Ir- 
relevant Clover (FLIC) fermion action in which the 
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FIG. 2: Mean square charge radius of the E + , as a func- 
tion of the pion mass squared, calculated in quenched lattice 
QCD 16]. The solid curve displays the finite volume fit, the 
infinite volume limit of quenched QCD is shown by the dashed 
curve, the valence sector of full QCD is represented by the 
long-dash-dot curve and the short-dash-dot curve shows the 
total result, including the disconnected loops. 
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FIG. 3: Absolute value of the mean square charge radius of 
the E~, as a function of pion mass squared, calculated in 
quenched lattice QCD |lq|. Curve descriptions are the same 
as those in Fig. 



irrelevant operators, introduced to remove fermion dou- 
blers and lattice spacing artifacts, are constructed with 
APE smeared links 17]. Perturbative renormalizations 
are small for smeared links and the mean-field improved 
coefficients used here are sufficient to remove 0(a) errors 
from the lattice fermion action [l8| . 

The (D(a)-improved conserved vector current 0] is 
used. Nonperturbative improvement is achieved via the 
FLIC procedure, where the terms of the Noether cur- 
rent having their origin in the irrelevant operators of the 
fermion action are constructed with mean-field improved 
APE smeared links. The results presented here are ob- 
tained using established techniques jjjj from a sample 
of 400, 20 3 x 40 mean-field 0(a 2 )-improved Luscher- 
Weisz |2l| gauge field configurations, having a lattice 
spacing of 0.128 fm, determined by the Sommer scale 
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FIG. 4: Mean square charge radius of a single u quark (of 
unit charge) in the proton, calculated in quenched lattice 
QCD [lg. Curve descriptions are the same as those in Fig. [5] 
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FIG. 5: Mean square charge radius of a single u quark (of unit 
charge) in the E + , calculated in quenched lattice QCD [T^. 
Curve descriptions are the same as those in Fig. |5] 



ro = 0.50 fen. Charge radii are obtained from dipole 
fits to quark-sector electric form factors, calculated at 
Q 2 = 0.227 ± 0.002 GeV 2 . The dipole form is well sup- 
ported physically and in lattice QCD simulations over a 
range of quark masses |22j. 

From Figs. [3 and |21 we see that the radius of the £ + 
grows significantly larger than that of the E _ as we ap- 
proach the chiral limit. This has a relatively simple, phys- 
ical interpretation in terms of the more extended spa- 
tial distribution of light quarks compared with strange. 
These simulation results arc fit using finite-range regu- 
larised (FRR), quenched chiral effective field theory for- 
mulated on a finite volume. Upon extracting the fit, one 
can restore the infinite-volume limit within the effective 
field theory — where corrections are observed to be quite 
small right down to the lightest simulated pion mass. 
The QCD valence sector and total QCD (including loop 
insertions) contributions are then estimated by replac- 
ing the long-range "tail" of QxPT by the corresponding 
xPT "tail" in the same fashion as Refs. 0, — which 



extends upon the empirical success observed in Ref. [23j . 

Using the techniques described in Ref. |f| , we estimate 
the effects of the dynamical sea and extrapolate the re- 
sults to the physical quark mass. The ratio of the S 
charge radii is found to be |E+/£ _ | = 1.141 ± 0.036 ± 
0.010 — where we observe the statistical fluctuations 
(first error) have been significantly reduced by construct- 
ing the ratio. The second uncertainty arises from the 
FRR scale dependence, A = 0.8 ± 0.1 GeV, constrained 
to provide agreement with experiment where available. 

In Figs. E| and we show the corresponding analysis 
for the u-quark mean square radii in the proton and S + , 
respectively. Here we observe that there is very little sen- 
sitivity to the mass of the spectator quark, with the ratio 
of the extrapolated valence contributions giving our best 
estimate of U p/u s = 0.993 ± 0.048 ± 0.000. The latter 
uncertainty indicates negligible sensitivity to the FRR 
scale. Putting these best estimates into Eq. we ob- 
tain e s = e R s d (l - fR^-^O.ll ± 0.06 ± 0.02 ± 0.33 fm 2 ], 
where the uncertainties are respectively, statistical, sys- 
tematic and experimental. 

In order to complete the evaluation, we also need to 
estimate the ratio of disconnected s to light quark loops, 
R s d - No attempt has yet been made to evaluate this ra- 
tio in an ab initio lattice QCD simulation, although fu- 
ture work could attempt this by building on the work of 
Refs. • We follow the same procedure as in Ref. 

but now for the charge radii rather than the magnetic 
moments, using the relative strength of the FRR loop 
contributions. With the FRR scale varying over the gen- 
erous range 0.8±0.2 GeV, this provides e R s d = 0.16±0.04. 

The final result for the proton strange, mean square 
charge radius is then (r 2 )f = ~ e s/3 = -0.007 ± 0.004 ± 
0.002 ± 0.021 fm 2 , where the uncertainties are statisti- 
cal, FRR scale and experimental in origin. Lattice scale 
determination uncertainties are negligible. At Q 2 = 
0.1 GeV 2 , we expect the first term in the expansion of the 
electric form factor to dominate, Ge(Q 2 ) = —Q 2 (r 2 )/Q. 
The strangeness electric form factor, conventionally de- 
fined without the charge factor, is then found to be 



G S E (0.1 GeV 2 



-0.009 ±0.005 ±0.003 ±0.027. (5) 



This is consistent with the analysis of world data at 
0.1 GeV 2 , namely G|;(0.1GcV 2 ) = -0.013 ± 0.028. The 
error on this theory calculation is of similar precision to 
the experiment. 

It is interesting that a combination of data from the 
SELEX collaboration at Fermilab, together with the con- 
straints of charge symmetry and modern lattice QCD can 
provide a tight constraint on the strange electric form 
factor of the proton. For the present time the errors are 
dominated by the experimental errors in the knowledge 
of the E~ charge radius. It would clearly be valuable to 
have new and more accurate measurements on both the 
£ _ and £ + hyperons. Data on the S~ (and S°) would, as 
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TABLE I: Results (fm 2 ) are compared with experiment. Un- 
certainties are statistical and systematic respectively. Note 
that (r 2 )f is reported including the strange-quark charge. 



Quantity 


This Work 


Experiment 


(r 2 ) P 


+0.685 ±0.047 ±0.051 


±0.757 ± 0.014 


(r 2 r 


-0.158 ±0.029 ±0.016 


-0.116 ±0.002 




-0.657 ± 0.037 ± 0.045 


-0.610 ±0.150 


(r 2 ) P s 


±0.001 ± 0.004 ± 0.004 


-0.010 ±0.022 



in the case of the strangeness magnetic moment, provide 
a valuable additional constraint. 

However, noting that the uncertainty in our result is 
dominated by the lack of precision in the observed S _ 
charge radius, we also consider a somewhat more ambi- 
tious approach. That is, we remove the dependence on 
the hyperons altogether and require only the observed 
nucleon charge radii and our best estimate for the quark- 
sector connected-current insertion, either u p or d p , 

In this case, uncertainties in the lattice scale determina- 
tion, a — 0.128(6) fm, become significant. Our best 
extrapolation of the valence sector contributions to the 
mean square charge radii in full QCD gives u p = 1.396 ± 
0.090±0.046±0.086 fm 2 and dP = 0.553±0.063±0.010± 
0.047 fm 2 , where the uncertainties are statistical, FRR 
scale and lattice scale in origin. These values provide two 
independent evaluations of the strangeness mean-square 
charge radius, (r 2 }§ = 0.000 ± 0.006 ± 0.007 fm 2 and 
(r 2 }f = 0.002 ± 0.004 ± 0.004 fm 2 , where the first uncer- 
tainty is the combined statistical and experimental er- 
rors, and the second is the combined FRR and lattice 
scale uncertainty. Combining these two in the extraction 
of the form factor at Q 2 — 0.1 GeV 2 gives 

G S E (0.1 GeV 2 ) = ±0.001 ± 0.004 ± 0.004 . (8) 

Following the procedure outlined in Ref. [2(|, charge- 
symmetry violations (CSV) of the maximum size esti- 
mated theoretically [2?J are found to be small with re- 
spect to the error on the experimental proton mean- 
square radius. The small magnitude of the strange form 
factors, predicted here and in Ref. || , suggests the exper- 
imental signal may be similar in size to the CSV effects. 
As the precision in experimental programs are improved 
a more careful QCD-analysis of CSV effects will be re- 
quired. 

Table H] provides a comparison of the present calcula- 
tion and experiment. The quoted errors on the present 
prediction are comparable to, or even better than, those 
anticipated in the forthcoming HAPPEx results. 
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